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We study numerically the out-of-equilibrium dynamics of interfaces at finite temperatures when driven well 
below the zero-temperature depinning threshold. We go further than previous analysis by including the most 
relevant non-equilibrium correction to the elastic Hamiltonian. We find that the relaxation dynamics towards 
the steady-state shows glassy behavior, aging and violation of the fluctuation-dissipation theorem. The interface 
roughness exponent a ~ 0.7 is found to be robust to temperature changes. We also study the instantaneous 
velocity signal in the low temperature regime and find long-range temporal correlations. We argue l//-noise 
arises from the merging of local thermally-activated avalanches of depinning events. 
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The dynamics of interfaces in random media has received 
much attention, both theoretical and experimental, in the last 
decade. This is largely due to the interest in driven inter- 
faces as models for nonlinear cooperative transport in many 
contexts, such as stochastic surface growth and kinetic rough- 
ening 1 1], magnetic flux lines in type-II superconductors |2|], 
the motion of charge density waves 01, propagation of frac- 
ture cracks 0], fluid imbibition in porous media 0], etc. The 
competition between elastic interactions and quenched disor- 
der leads to a rich behavior, which is not yet fully understood. 
A key macroscopic observable to look at is the average inter- 
face velocity v as a function of the external driving force /. At 
zero temperature, quenched disorder tends to pin the interface 
and this leads to a depinning transition from a pinned phase 
(v = 0) to a moving phase (v > 0) at a critical value of the 
external driving force f c . However, one particularly important 
question is how the system responds to this external forcing in 
the presence of non-negligible thermal fluctuations. 

At finite temperatures, when the interface is driven well be- 
low the zero-temperature critical force f c (T = 0), the inter- 
face velocity v is extremely small. The interface then moves 
in bursts of activity, which are thermally activated and, in the 
steady state, the response of the system to a small external 
force exhibits glassy properties known as creeping 0]. Ex- 
periments 0] and theoretical studies 0] have confirmed the 
existence the creep behavior. Theoretical approaches to this 
problem using mean-field or functional renormalization group 
techniques 
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are usually valid in high di- 
mensions and difficult to extend to one dimensional interfaces. 
More importantly, these approaches treat the dynamics of an 
elastic string in a random potential, i.e. interaction terms de- 
rive from an elastic energy H = (v/2) J dx[*J (V/i) 2 + 1 + 
$ p (.t, h)] in a random pinning potential <P p in the small slopes 
limit |V/i| -C 1. However, the elastic line approach gives 
a roughness exponent a m 1.2 at the depinning transition, 
which leads to difficulties in the thermodynamic limit since 
interface fluctuations ([h(x) — h(x)] 2 ) ~ L 2a diverge and, 
indeed, the small slopes approximation breaks down Il4jll51l . 
Recent theoretical arguments have pointed out that a 



Kardar-Parisi-Zhang (KPZ) nonlinearity flrjll may emerge 
from anisotropy of the disorder 11711 and/or from higher or- 
der terms in an elastic expansion Il4l fl^l . Also, recent de- 
velopments in imaging technology have allowed to directly 
visualize the boundary between the vortex invaded and the 
vortex free (Meissner) regions in superconductors. In a se- 
ries of nice experiments Surdeanu et al. fUjl recently mea- 
sured the temporal and spatial correlations of these fronts for 
a YBa2Cu307_a; superconductor. They reported interfaces 
that exhibited scale-invariant kinetic roughening with expo- 
nents consistent with KPZ dynamics. Similar results were 
later obtained in Ref. fl^l for thin Nb superconductor films. 
These theoretical and experimental results make evident the 
importance of understanding the effect of thermal fluctuations 
on pinned interfaces when the interaction between degrees of 
freedom contains nonequilibrium contributions, like the KPZ 
nonlinearity. 

In this Letter we study the relaxation towards the steady 
state of an interface locally pinned by a disordered back- 
ground. Our study includes the most relevant KPZ non- 
equilibrium correction to the elastic energy. We focuss on the 
interface dynamics well below the depinning transition, when 
only temperature can locally detach the interface. We find that 
the system very slowly relaxes (creeping) to a steady-state. 
During the relaxation regime the system exhibits aging and vi- 
olation of the fluctuation-dissipation theorem. The roughness 
exponent of the string is found to be a w 0.7 at short scales 
and robust to changes in temperature. Remarkably, the instan- 
taneous velocity signal exhibits long-range temporal correla- 
tions. We argue l//-noise arises from the merging of local 
thermally-activated avalanches of depinning events. 

We consider a ID line described by a single-valued func- 
tion h(x, t), giving the transversal position h of the front from 
the h = axis at time t. The front is moving through a het- 
erogeneous medium, which can be described by a quenched 
disorder rj(x, h) and the interface obeys the equation of mo- 
tion: 

d t h = v V 2 h + A (V/i) 2 + VT £(x, t) + a T](x, h) + f, (1) 
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FIG. 1 : (Color online) Average height as a function of time for the 
nonlinear equation. The different curves correspond from bottom to 
top to the temperatures T = 0,0.01,0.04,0.16 and 0.64. In the 
inset, the spatial power spectra of the interfaces for the three lowest 
temperatures (colors match in both graphs). The dashed lines have 
slope 1 in the main plot and slope — (2a + 1) = —2.4 in the inset. 



where t) is a gaussian white noise with zero mean and 
unit variance describing the thermal fluctuations at temper- 
ature T, v is the elastic constant, and the KPZ nonlinear- 
ity A (Vft) 2 gives the most relevant nonequilibrium correc- 
tion to the equilibrium elastic energy. We only consider 
here the case of non-correlated gaussian disorder so that 
(t}(x, y)v(x', y') = S(x - x') 6{y - y') and (r)(x, y)) = 0. 

Interface scaling.- We have carried out extensive numer- 
ical simulations of the Langevin equation of motion 0. To 
solve Q numerically we discretized h(x,t) along the x direc- 
tion with lattice spacing Ax = 1 and used a stochastic Eu- 
ler scheme to integrate the equation of motion with periodic 
boundary conditions. The constants v and a can be scaled out 
and without lose of generality we have set v = a = 1. All 
the results reported in the following were obtained for an in- 
tegration time step At = 0.01, which was enough to guaranty 
stability of the numerical scheme and Ai-independence of the 
numerical results. We are interested in the fully nonlinear 
regime, thus, for the system sizes we used (up to L = 2046), 
the strength of the nonlinear term is set to A = 1. Our results 
do not depend on A, as far as it is large enough to allow the 
nonlinearity to fully operate well before saturation is reached. 
For comparison we have also studied the usual elastic regime 
for A = 0. 

Simulations at T = revealed that the depining transi- 
tion takes place at f c (T = 0, A = 1) = 0.36 ± 0.01 and 
f c (T = 0,A = 0) = 0.99 ± 0.01 for the nonlinear and 
linear models, respectively. A more precise determination of 
the critical thresholds is not required since we are interested 
in the dynamics of the system well below the critical point. 
In the following, unless otherwise stated, we will always be 
driving the system with external forces /(A = 1) = and 
/(A = 0) = 0.1, well below the critical threshold. We now 
focuss on the effects of varying the external temperature T 



when the system is in the zero-temperature pinned phase. 

We always start our simulation with a flat interface h = 
and then let the system evolve at a given temperature T. 
We monitor the average height (h) from which the inter- 
face velocity can be obtained in the long times limit v = 
]xm t ->oo{h) /t. 

Initially, as may be seen in Figure 1, the front advances 
very fast. The interface moves through the quenched disor- 
der and progressively gets locally trapped at more and more 
pinning sites. The duration of this transient is system size in- 
dependent and relatively short. After that initial outburst, the 
interface begins relaxation towards the steady state. The latter 
is characterized by a constant velocity. The time span of the 
relaxation process strongly depends on temperature, increas- 
ing very dramatically at low T. During relaxation, the average 
height displays a nonlinear growth in time, (ft) ~ t 7 ^ T \ char- 
acteristic of the creeping regime. The actual value of 7(T) 
decays rapidly with decreasing T, being the average height at 
very low T compatible with a logarithmic growth. 

This singular dynamics of the average height during the 
creeping regime imposes constrains upon the evolution of the 
interface width. Since w(t, L) ~ t a / z cannot grow faster than 
the height (ft) ~ f ', we have ct/z < j(T). Simulations show 
that the value of the roughness exponent a w 0.7 remains 
essentially constant at low temperatures. In the inset of Fig- 
ure 1 we plot the structure factor S(k, t) = (h(k, t)h(—k, t)) 
at long times t = 10 4 for three different temperatures T = 
0, 0.01 and 0.04 (where h is the the Fourier transform of the 
interface height in a system of lateral size L, with k being 
the spatial frequency in reciprocal space). The inset clearly 
shows that (i) after long times 10 4 the interface has not yet 
reached a steady state, (ii) the interface is scale-invariant 
S(k,t) ~ fc-( 2a+1 ) at short scales k > t^ 1/z , and (Hi) the 
roughness exponent a does not depend on temperature. This 
in turn implies that it is the dynamic exponent z the one chang- 
ing with temperature (since a/z < 7(T)). This is made more 
evident as T — > 0, when dynamics is dramatically slowed 
down and z suffers a significative increment. 

Aging.- Ultra-slow relaxation towards the steady state in- 
dicates typical out of equilibrium glassy features in the sys- 
tem, in particular the existence of aging-namely, the break- 
down of the time translation invariance symmetry I^ EolEHl . 
In order to study aging in our system we compute the two- 
times connected correlation function 

r(t, t w ) = ([h(x,t) - (h)(t)\ [h{x,t w )~ (ft) (*„)]). (2) 

If the system ages, T(t, t w ) must be a function of both t and 
t w , and not only of their difference t — t w . We follow the 
standard procedure, which closely mimics what can be done 
in real experiments. The interface is evolved at some very high 
fixed temperature (we typically used T = 4 for the nonlinear 
equation and T = 36 for the linear case) to anneal the system 
until the stationary state is reached. Then, when the system 
is in the steady-state, the temperature is suddenly decreased 
to the desired measuring level T. This sets up the system out 
of equilibrium and the time counter is started, t = 0, at the 
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FIG. 2: (Color online) Connected correlation as a function of the 
time difference t — t w . The main graph is for the nonlinear equation 
at T = 0.04 and t w = 250, 500, 1000 and 2000 from left to right. 
In the inset, the same plot but for the linear equation, at T = 0.49 
and the same set of t w values. 




FIG. 3: (Color online) Parametric representation of the susceptibility 
as a function of the connected correlations for the nonlinear model 
with T — 0.01 (black curves) and 0.02 (red ones). The waiting times 
in both cases are t w = 500 (circles) and t w = 1000 (squares). 



freezing event. We then monitor the relaxation dynamics by 
measuring T(t, t w ). Typical T(t, t w ) curves for a set of wait- 
ing times t w are shown in Figure 2 for a external temperature 
T = 0.04 and 0.49 for the nonlinear and linear cases, respec- 
tively. Both the nonlinear and the linear equations show clear 
aging during relaxation towards the steady state similarly to 
other glassy systems. 

Fluctuation-dissipation violation.- Another important as- 
pect of glassy dynamics is the way in which the fluctuation- 
dissipation relation is broken. At equilibrium, the suscepti- 
bility is linearly related to the spatial correlations of the con- 
jugated variable. A modification for slowly evolving out of 
equilibrium systems has been recently proposed by Cuglian- 
dolo, Kurchan and Peliti lilll . through the generalized re- 
sponse function 

X[T] dT(t,t w ) 
R(t,t w )-— , (3) 



where R(t,t w ) = (Sh(x,t)/5f(x,t w )), i.e., the response of 
the conjugated variable (the interface height in our case) at a 
time t to the application of a local external field (the driving 
force / in our case) at the same position x at a previous time 
t w . X[T] is a functional that can vary between two asymptotic 
constant values, either it takes value one (for scales that are 
already in the stationary regime) or X[T] — > X^ (for those 
scales that are still out of equilibrium). 

Numerically, it is more convenient to study the suscepti- 
bility, defined as x{t;tw) = f+ dsR(t,s), instead of the 
response function R(t, t w ). The reason being that the com- 
putation of x(~b>~kw) does not require the introduction of an 
instantaneous spike-like perturbation in the external field and, 
in addition, the integration in time implies that the external 
field operates during a period of time making the estimation 
of its effects much easier to detect. We integrate R(t, s) in Eq. 



l|3} in the variable s and obtain 

X(t,t w ) = ™\r(t,t)-T(t,t w )]. (4) 

In order to check whether the fluctuation-dissipation relation 
is broken during relaxation in our system, we use a similar 
numerical setup as before. We anneal the system by heating 
at very high temperatures until a steady-state is reached, then 
the temperature is decreased and that event marks the time 
origin t = 0. After a waiting time t w , we run in parallel 
two copies A and B of the interface. Copy A is then driven 
with constant force f A (x,t) = f A as before, and copy B is 
driven with a slightly different local driving force f B (x, t) = 
f A + Sf(x), where 8f(x) is a small uncorrected random field 
which takes at random the values ±e with equal probability. 
The susceptibility can then be approximated by x(t>tw) ~ 
((h B (x,t) — h A (x, t)) Sf(x))/\e\. In Figure 3, parameteric 
plots with the susceptibility as a function of T(t, t) — T(t, t w ) 
are shown for the nonlinear equation (A = 1) and for different 
temperatures. The fluctuation-dissipation theorem is clearly 
violated, again akin to observations in other glassy systems. 

Thermally activated avalanches.- Finally, we pay atten- 
tion to the dynamics of avalanches at finite temperatures in 
the nonlinear model Eq. Q. Avalanches in forced supercon- 
ductors have attracted much attention I23I E4I1 . but little is 
known about avalanche dynamics at a finite T. h(t) data (not 
shown) for a single disorder realization have a dented and 
devil staircase-like appearance caused by the stick-slip mo- 
tion associated with avalanches triggered by a local thermally- 
activated depinning event. We monitor the avalanche dynam- 
ics in our system by looking at the temporal correlations of the 
average velocity signal v(t) = dth(t) well below f c (T = 0). 
In Figure 4, we plot the spectral density S(w) = (\v(w)\ 2 ), 
averaged over disorder realizations. The velocity signal v(t) 
shows long-range temporal correlations. In the low-frequency 
limit we observe a power-law decay S(w) ~ w~ e with a 
temperature-dependent exponent 9(T) ranging from uncorre- 



4 




log w Ceo) 



tures (T ss 0.25 in Figure 4) all the pinning sites are overcome 
by the merging of local avalanches that can extend across the 
whole system freeing the interface from the disorder. Further 
increasing the temperature above that point will not affect the 
dynamics since temperature is high enough to overcome local 
pinning and the depinning transition is actually wiped out. In 
this regime the interface is described by the KPZ dynamics 
and w -1 / 3 velocity-fluctuation spectra. 
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No. 0312510 (USA) and by the Ministerio de Education y 
Ciencia (Spain) under projects BFM2003-07749-C05-03 and 
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FIG. 4: (Color online) Spectral densities of the velocity signal for 
the nonlinear equation (A = 1). Temperatures are T = 0.01, 0.09, 
0.1225, 0.16, 0.25, 1 and 2.25 from bottom to top. The straight lines 
are guide for the eye and have a slope —0.3 and —1.16 respectively. 

lated noise 9 ~ (at very low temperatures) to long-range 
correlations 1.0 < 9(T) < 1.2 for temperatures in the range 
0.05 < T < 0.20. For even higher temperatures, T > 0.25, 
velocity spectra clearly show a crossover to a different power- 
law decay with exponent 6(T) « 0.30, which then remains 
unaltered when temperature is further increased. The latter 
regime corresponds to a free-moving interface in the KPZ uni- 
versality class, which velocity-fluctuation spectra is known to 
diverge as 1 jw e at low frequencies, where = (d + 4)/z — 3 
in d + 1 dimensions I25I1 . In d = 1 one then finds that 6 = 1/3 
in good agreement with our numerical results at high temper- 
atures where the interface is freely moving with a finite dis- 
placement velocity. We obtained similar results for the linear 
case (A = 0), where 6 ~ 1 in the intermediate temperatures 
regime, and the crossover to w° for the highest temperatures, 
as expected for the linear equation I25I1 . 

This change in the temporal correlations of the velocity sig- 
nal at high enough temperatures can be pinned down to the ex- 
istence of thermally activated avalanches of depinning events. 
When the interface is driven well below f c at finite tempera- 
tures, uncorrected avalanches of typical size £ are triggered 
all over the front. At very low temperatures thermally acti- 
vated avalanches are very unlikely and localized events, dis- 
joint and uncorrected to one another, giving rise to a typi- 
cal flat spectrum (8 ~ 0) of the velocity signal. However, 
at higher temperatures we expect more sites to detach from 
disorder and so more local avalanches to occur. Avalanches 
can overlap and merge to give rise to a larger avalanche, even- 
tually covering a macroscopic region of the interface. This 
leads to bursts of spatially connected moving events and, con- 
sequently, long-range temporal correlations (6 ~ 1) in the 
interface velocity fluctuations. This can be seen in Figure 4 
for 0.09 < T < 0.16. The higher the temperature, the more 
likely a site gets depinned, triggering a local avalanche of de- 
pinning events of size £. Eventually, at high enough tempera- 
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